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SUMMARY 

Local turbulent heat -transfer coefficients were obtained experimentally for 
the dissociating system N 2 0 4 -N0 2 flowing through an electrically heated tube. 

For these experiments the tube length to diameter ratio was 40. Experimental 
conditions covered a seven-fold increase in the uniform heat flux up to 44 000 
Btu per hour per square foot and an eight-fold increase in mass average velocity 
up to 130 000 pounds per hour per square foot. The maximum tube wall tempera- 
tures ranged from 359° to 572° K. 

Modified local heat-transfer coefficients, based on wall- to bulk-enthalpy 
differences, were correlated very well as a simple power of the mass average 
velocity. 

The local modified heat -transfer coefficients were correlated with a Dittus- 
Boelter -Colburn-type relation with all equilibrium transport and thermal 
properties evaluated on a mean temperature integral basis. The best curve 
through the data was about ten percent below the curve from constant property- 
analyses . 


INTRODUCTION 

The absorption of the heat of chemical reaction by a dissociating gas is 
somewhat analogous to the absorption of the heat of vaporization of a boiling 
liquid in that both fluids act as effective heat sinks. For this reason a num- 
ber of gases have recently been examined as potential heat-transfer fluids, 
particularly for special applications in cooling of nuclear reactors (ref. l) . 

In a heat-transfer process the effect of a chemical reaction is most appar- 
ent in the unusual sensitivity of the fluid thermal conductivity and specific 
heat to changes in temperature and pressure. Calculation of these properties is 
very difficult except when the products and reactants are in chemical equilib- 
rium. For this special condition the expressions for the thermal conductivity 
and specific heat are relatively simple and tractable (ref. 2). 



Of the possible equilibrium chemical reactions, the gaseous system 
NgOq 2N0 2 has been most attractive for heat-transfer studies because its in- 
teresting thermophysical properties can be exploited in experiments -within the 
temperature range 300° to 500° K and at pressures of 1 to 2 atmospheres. The 
gas can also be obtained and stored in suitably pure form. The Ng 0 4 system was 
used in the present experiment because of both the previous reasons and the re- 
cent reported study of its transport properties (ref. 3), which gave necessary 
data for the correlations . 

The large variations of thermal properties in a reacting system provide a 
severe test for the standard turbulent flow heat-transfer correlations, which 
have usually been derived for constant property fluids. The purpose of the 
present experiment was to provide extensive heat-transfer data for the turbu- 
lent flow in a tube of a gas in chemical equilibrium and with these data check 
some of the correlations that have been devised to account for deviations from 
the constant property assumptions in the analyses. 

In the present experiments the heat exchange was from an electrically heated 
tube to turbulent flows of gas. The experimental conditions covered a 

seven-fold increase in the uniform heat flux up to 44 000 Btu per hour per 
square foot and an eight-fold increase in the mass average velocity up to 
130 000 pounds per hour per square foot. The tube wall temperatures were ad- 
justed to provide wall properties with maximum chemical reaction effects for 
the low wall temperature runs and wall properties due to the completely dis- 
sociated gas for the highest wall temperature runs. The experiments were de- 
signed to give the largest variations in the important heat-transfer variables 
and parameters as was possible with the equipment. 

Sources of data from previous investigations of the turbulent heat- 

transfer system are listed in table I. References 4 to 7 reported ex- 
perimental results with uniform wall temperature conditions. These experiments 
are characterized by relatively modest values of wall heat flux. 

The other references in table I are concerned with heat transfer under uni- 
form heat flux conditions. The data of reference 8 were taken at the very short 
length to diameter ratio of 2 and were only later correlated (ref. 9) using a 
Colburn-type equation and more complete thermophysical property data. 

Two companion papers (refs. 10 and ll) considered radial variations in prop- 
erties in correlating the data with both Deissler's analogy (ref. 12) and with 
Colburn's analogy (ref. 13, p. 424). There was fair agreement between theory 
and experiment. 

More general discussions of reacting systems with heat transfer, and with 
emphasis on boundary layer flows, have been presented by Spalding (ref. 14) and 
Knuth (ref. 15). 
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SYMBOLS 


A f p cross-sectional flow area, sq ft 
inside heated wall area ; sq ft 

A 

Cp mass basis specific heat, Btu/(lb)(°F) 

Cp^ e mass basis equilibrium specific heat, Btu/(lb)(°F) 

0p^ e mean temperature integral specific heat, Btu/(lb)(°F) 

Cp^-p mass basis frozen specific heat, Btu/(lb)(°F) 

✓S. 

Cp^ r mass basis reactive specific heat defined by eq. (A27) 

D inside tube diameter, ft 

@ binary diffusion coefficient, sq ft/hr 

E electrical potential, V 

AF° standard state Gibbs 1 free energy, cal/g mole 

f friction factor 

G mass velocity, W/A^, lb/(sq ft) (hr) 

H° standard state molar enthalpy, cal/g mole 

AH° standard state heat of reaction, cal/g mole 

AH° standard state heat of formation, cal/g mole 

h heat-transfer coefficient based on wall to bulk temperature difference, 

Btu/ (hr) (sq ft)(°F) 

h 1 modified heat -transfer coefficient based on wall to bulk enthalpy dif- 

ference, lb/(hr)(sq ft) 

I electrical current, A 

i mass basis enthalpy, Btu/lb 

Colburn factor, St T (Pr e ) 2 / 3 

Kp equilibrium constant for system, atm 
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equilibrium thermal conductivity, Btu/(hr)(sq ft)(°F/ft) 
mean temperature integral thermal conductivity, Btu/(hr)(sq ft)(°F/ft) 
frozen property thermal conductivity, Btu/(hr)(sq ft)(°F/ft) 
reactive contribution to thermal conductivity, eq. (A30) 
heated tube length, ft 

molecular weight of equilibrium gas mixture, Ib/lb mole 
molecular weight of NOg molecule, 46.008 lb/lb mole 

molecular weight of NgO^ molecule, 92.016 Ib/lb mole 

mass fraction of species 
coefficient in eq. (7) 

standard Nusselt number, hD/k^ or hD/k e 
modified Nusselt number, h*D/r e 

modified Nusselt number based on mean temperature integral equilibrium 
properties, h*D/P e 

number of fillings of test loop in eq. (l) 
absolute total pressure, atm 
Prandtl number, p/r 

Prandtl nuunber based on equilibrium properties, p/r„ = pC /k 

Prandtl number based on mean temperature integral equilibrium proper- 
ties, p/r\T 

Prandtl number based on frozen properties , p/r f 

absolute partial pressure, atm 
rate of electrical heat input, Btu/hr 
heat flux to gas, Btu/(hr)(sq ft) 



R universal gas constant 

Re Reynolds number, GD/j jl 

Re Reynolds number based on mean temperature integral viscosity, GD/jT 
Bt T modified Stanton number, Nu*/R e Pr e or h*/G 
T temperature, °K 

U bulk velocity, G/p^, ft/hr 

W mass flow rate, Ib/hr 

x distance from entrance of heated test section, ft 
y mol fraction, mole/mole mixture 

a degree of dissociation of initial 1 mol ^ole/mole 

a 

r thermal exchange coefficient, k/C^ 

T equilibrium thermal exchange coefficient, k^/fi , lb/(ft)(hr) 

e c p , e 

r e mea n tem perature integral equilibrium thermal exchange coefficient, 

p ,lb/(ft)(hr) 

P frozen thermal exchange coefficient, kp/C ,lb/(ft)(hr) 
i 1 P j i 

p dynamic viscosity, lb/(ft)(hr) 

p mean temperature integral dynamic viscosity, lb/(ft)(hr) 

p mixture density, Ib/cu ft 

cp general transport property in eq. (17) 

X Boltzmann factor, 1.4384 cd/t 

u> wave number, cm^ 

Subscripts : 
b bulk 

e equilibrium 

f frozen, or constant composition 


mix 


mixture 



P refers to properties at constant pressure 

ref reference (indicates evaluation of properties at reference temperature 

T or at temperature corresponding to reference enthalpy i ) 
rex ref 

w wall 

0 refers to property at absolute zero temperature 

Superscript: 

mean temperature integral property 


APPARATUS AND PROCEDURE 
Flow System 

The apparatus formed a closed flow loop containing a compressor, two ro- 
tameters, an electrically heated test section, and water-cooled heat exchangers 
as basic elements. The system is shown schematically in figure 1. The entire 
rig, exclusive of the control panels, was enclosed within a fume hood fitted 
with Lucite panels. 

The compressor had a rated capacity of 17.0 standard cubic feet per minute 
at a maximum pressure rise of 15 pounds per square inch gage. 

Two rotameters connected in parallel provided the gas flow rate indication. 
With both rotameters fully open the maximum indicated flow rate was 279 pounds 
per hour of EgO^ gas at a standard density of 0.385 pound per cubic foot. 

Two sets of jacketed heat exchangers, one upstream and the other downstream 
of the compressor, were fitted with both cold and hot water inputs so that the 
cooling water temperature could always be kept above the dew point of the test 
gas in the loop to avoid condensation. 


Power Supply 

Electrical power to the test section was supplied by a regulated saturable 
core reactor with a rated capacity of 25 kilovolt amperes and a maximum output 
of 10 volts at 60 cycles on the transformer secondary. 


Test Section 

Details of the test section and mixing boxes are given in figure 2. The 
test section was a tube of 347 stainless steel, with a 0.500-inch outside diam- 
eter, 0.404-inch inside diameter, and an overall length of 20 inches between 
mixing box flanges. The first 4 inches of the test section were unheated. Two 
pressure taps were provided on the test section, one on the upstream electrical 
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connection, and the other on the exit flange. The mixing boxes were constructed 
basically as three concentric cylindrical cans, which provided four abrupt turns 
with the mixing of the fluid stream. The mixing tanks were fabricated from 
347 stainless -steel sheet. 

3 

Although not shown in figures 1 and 2, a 20-inch-long transite tube, Sc- 
inch inside diameter, was fitted around the test section between the entrance 
and exit flanges. The annular space between the test section and the transite 
tube was filled with loose insulation. The transite tube was cut in half longi- 
tudinally so that it could be fit around the test section. Test section thermo- 
couple leads were also brought out through the longitudinal joint. 


Ins tr ument at i on 

Temperatures were sensed with 24-gage iron -const ant an thermocouples and 
were recorded on a flight recorder. Thermocouples were swaged in sealed 
stainless -steel sheaths where direct contact with the gas was necessary (e.g., 
at the two rotameters), in the inlet and outlet of the compressor, and in the 
mixing tank gas streams at the inlet and outlet to the test section. 

Bare thermocouples were spot welded to the outside of the test section, as 
is indicated in figure 2. These were usually put on in pairs, on the top and 
bottom, at any one longitudinal position. Each individual wire was looped once 
around the test section, and was held off the tube surface with insulating 
beads, before being brought out to a junction box. Five thermocouples in the 
middle 6 inches of the heated test section were alternated top and bottom. In 
total there were 30 thermocouples on the test section, and 22 on the transite 
insulating shell. The thermocouples on the transite shells were evenly distri- 
buted along its length, alternating top and bottom. 

In addition to the flight recorder, thermocouple temperatures were indica- 
ted on a precision potentiometer in conjunction with an electronic null point 
indicator. The potentiometer temperatures were taken relative to the ice point. 
However, it was also possible to take temperature differences directly with the 
potentiometer. These more precise readings were used only for the crucial mix- 
ing tank and rotameter temperatures. 

It is indicated in figure 1 that static pressures at the rotameter inlets 
and at the test section were obtained on manometers. Pressures at the inlet and 
outlet of the compressor were not needed with much precision and thus were mon- 
itored by stainless steel Bourdon gages. 

All static pressure tubing was of l/8-inch thinwall 347 stainless steel and 
was always attached to the top of the stainless drip pots. These pots were nec- 
essary to trap any gas from the system that may have condensed at the ambient 
temperatures. This ensured that the manometer lines were always free of liquid. 

The gravity pots separating the manometers from the drip pots, which are 
shown in the lower right-hand corner of figure 1, were a later addition to the 
equipment to deal with the problem of HOg gas entering and reacting with the 
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manometer fluids. These gravity pots, tvo in number in each manometer line, 

were welded stainless -steel cans about 6 ^ inches high and 2^- inches in diameter. 

Each set of pots in series functioned as a trap for the EOg gas because the gas 

had a density of about three times that of air at the ambient conditions of the 
experimental runs. 

Pressure drops across the test section were obtained with a differential 
manometer . 

The voltage drops across the test section were measured directly on lab- 
oratory precision voltmeters. Current was measured by a conventional tech- 
nique: an ammeter in conjunction with a current transformer having a ratio of 

160. A wattmeter was also used with the voltmeter and ammeter as a check on 
their accuracy. 


Corrosion Problems 

Special precautions were necessary to counter corrosion by gaseous NOg. 

The flow loop was constructed as far as it was possible of glass or stainless 
steel. Gaskets were constructed from sheet Teflon, and all screwed fittings 
were doped with a Teflon-water paste; rotameter tube washers were machined from 
Teflon rod; compressor shaft seals and bearings were lubricated with Kel-F oil. 

The cast iron body of the compressor and the plastic vanes of the rotor 
were not attacked by the gas. However, the compressor had to be operated dry 
(i.e., without any vane lubricant) after a number of oils apparently were po- 
lymerized and degraded by the gas. 

All bare iron-constantan thermocouple junctions corroded on first con- 
tact with the gas. Thereafter, it was necessary that all thermocouple probes 
in direct contact with the gas be constructed with the thermocouple wires swaged 
in sealed stainless -steel sheaths. 


Filling and Operation 

Filling the test loop with the N0 2 gas was accompanied by the purging of 
the air in the system. This was accomplished by the pulling of a modest vacuum 
on the system followed by filling of the loop to some specified pressure from a 
small pressurized supply bottle of N0 2 exterior to the system. The supply bot- 
tle was not stored in the vicinity of the rig, and is not shown in figure 1. 
Running the compressor for 10 minutes after a system filling with N0 2 was con- 
sidered sufficient to thoroughly mix the N0 2 gas with residual air in the loop. 
Each cycle of evacuation, filling, and mixing reduced the mole fraction of 
residual air in the system. 

Since the equipment was always filled to the same pressure and evacuated 
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to the same pressure, then the mole fraction of residual air after n cycles 
was 

( 1 ) 



That the air was rapidly exhausted from the system can he seen from a nu- 
merical example. The loop was typically exhausted to 25 inches of vacuum (i.e., 
p = 5/30 =l/6 atm abs) and filled to 2 atmospheres absolute. Then after 
3 cycles 





0.00058 


This number represents a nominal impurity of 0.058 mole percent after 3 cycles, 
which is well within the 99.5 weight percent minimum assay claimed for the cyl- 
inder gas (ref. 16) . 


After the filling operation, power was applied to the test section for a 
given gas flow rate and the test section with its blanket of insulation was al- 
lowed to come to thermal equilibrium, which usually required 4 to 5 hours. The 
compressor shaft seals never worked perfectly and allowed a small but steady 
leak of N0 2 gas into the fume hood. This meant that an equivalent amount of 
Wg0 4 was constantly fed to the system in order to maintain constant pressure. 

The water temperature in the cooling jackets of the heat exchangers was 
adjusted according to the vapor pressure curve given in figure 3. For a given 
system pressure, the temperatures at all points in the loop had to be at values 
above the vapor pressure curve. This put a lower limit on the amount of cooling 
of gas from the hot test section, and a limit on the amount of recooling of the 
gas heated in the compressor. The penalty for having liquid N2O4 in the loop 
was flow and temperature surges . 

A chemical analysis or assay of the oxides of nitrogen in the system was 
made during the time that the rig was coming to thermal equilibrium prior to a 
run. The complete procedure for the assay is given in reference 16. Only a 
resume will be sketched here. 

Three clean, dry glass balloons (e.g., Guillard bulbs) of about 500- 
milliliter capacity were evacuated and weighed on an analytical balance. The 
bulbs were filled with gas from the test loop to about 1-atmosphere pressure and 
then weighed again. A measured amout of a dilute H 2 0 2 solution was forced into 

the bulbs. This oxidized any lower nitrogen oxides to W0 2 , and this gas dis- 
solved in the solution to give HNO3. The bulbs were opened, and the acid was 
decanted and washed into Erlenmeyer flasks where it was titrated to a methyl red 
endpoint with standard 0.2N NaOH reagent. The amount of KNO3 determined in each 
Guillard bulb, expressed on a weight basis as W0 2 , was related on a percent ba- 
sis to the total weight of original gas sample. The accuracy of the analysis 
was not exceptional, about 0.2 percent. However, the assay indicated that air 
was essentially exhausted from the system, for the weight percent of N0 2 was 

never lower than 99.5, which was the minimum manufacturer's specification. 
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DATA REDUCTION AND CORRELATION 


The raw data from the experiments consisted of the voltage and current 
readings; the indicated temperatures at the rotameters, mixing boxes, test sec- 
tion walls, and transite insulation walls; the barometric pressure; the manom- 
eter pressures for the rotameters and the ends of the test section; and the in- 
dicated rotameter readings. The total power dissipated in the test section was 


Q = (3.143) (160)EI (2) 

where the voltage E and current I were values corrected with meter calibra- 
tion curves. 

Overall heat loss was obtained by comparing the integral mean wall tempera- 
ture with the corresponding mean temperature from the static heat loss tests. 
This provided a corrected heat input, 

^corr = ^ %oss 

This value was compared with the total enthalpy rise of the test gas between 
the mixing boxes to obtain a check on the heat balance, 

tt . -1 / ,\ ^corr " ^^b, exit ’ ^inlet^ __ /„\ 

Heat balance error (percent] =* ^ — — X 100 (3] 

^corr 

Temperatures on the inside of the wall were calculated from the external 
thermocouple readings using equation (l) of reference 17. The maximum tempera- 
ture drop through the wall was 3° K, which occurred at the maximum heat input of 
6240 Btu per hour. A summary of the important experimental test data is given 
in table II. 

From equation (2) of reference 17, the heat input to the test gas was com- 
puted from 

Q ~ -^b, exit “ ^b, inlet) 

and these are the values which have been given in table II. The bulk enthalpy- 
rise values of heat input were used in subsequent calculations. The measured 
electrical heat inputs were subjected to nonsinusoidal variations in waveform 
of the power supply. These waveform variations did not affect the bulk temper- 
ature rise, but they did give trouble in the recording of the electrical meter 
readings. The erratic waveform was the principal reason for the variations in 
the heat balance error recorded in table II. 

The double thermocouple probes in each mixing box agreed to within 0.1° R 
for all tests, which was also the estimated accuracy of measurement of these 
temperatures . 

The test rig was given a preliminary checkout of procedure and accuracy 
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with air as the test gas* The excellent results obtained, which proved the 
equipment, have been briefly reported elsewhere (ref. 18). 


PROPERTIES 

The methods used to compute the thermodynamic and transport properties of 
the gaseous NgO^ system are discussed in the appendix. 

The magnitude of the effect of the chemical reaction on the thermal con- 
ductivity and specific heat of the gas mixture can be seen from the analytical 
values in figures 4 and 5. The difference between the "frozen" (i.e., constant 
composition) and equilibrium curves is the measure of the effect of reaction. 
The effect diminishes rapidly at higher temperatures where the system is essen- 
tially dissociated to NOg. 

The absolute enthalpy of the system, shown in figure 6, exhibits an in- 
flection point at about 340° K for the case of P * 2 atmospheres. As can be 
deduced from the material in the appendix, this represents the neighborhood of 
the maxima for the equilibrium k e and c 

Pj e 

Dynamic viscosity of the system, shown in figure 7, is much less affected 
by the presence of chemical reaction. It is, however, composition dependent as 
can be seen by the two analytical curves for one and two atmospheres. In this 
regard viscosity behavior is very similar to that of frozen thermal conductiv- 
ity. In fact, the methods of computation of kp and q are similar. 

The thermal exchange coefficient T and Prandtl number Pr shown in fig- 
ures 8 and 9, respectively, are properties derived after their respective com- 
ponents have been computed. 


DATA CORRELATION 


The conventional correlation of fully developed turbulent flow heat trans- 
fer in a tube is the Dittus-Boelter equation (ref. 19). 


hD 

k 



(5) 


For a uniform heat flux to a nonreactive fluid, the fully developed heat- 
transfer condition has a constant heat-transfer coefficient 

( 6 ) 



which demands the T w - T-^ term be constant. 

Humble (ref. 17) and Weiland (ref. 20) have discussed the important prob- 
lem of choosing the temperature at which the transport properties are calcula- 
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ted for use in equation (5) . Wall and bulk temperatures are the most conven- 
ient to obtain, while the best results seem to come from the use of a reference 
temperature defined by 

( 7 ) 

with the value of m usually around 0.5 for nonreactive gases (ref. 13, p. 393). 



The uniform heat flux of equation (6) can be related to the mixing box 
temperatures by 


where 


Q qA w - ¥ C p ^ b } exit “ T b, inlet j 



dx 


( 8 ) 


( 9 ) 


Equation (8) is equivalent to the enthalpy expression 


^ W ^b,exit ~ i b,inlet^ 


( 10 ) 


provided that the correct specific heats are used in equation (9). For a chem- 
ical equilibrium situation like the present system, the use of frozen 

Cp^p in equation (9) can give values of predicted Q from equation (8) that 
are almost an order of magnitude lower than the actual measured heat input. 


Since equilibrium enthalpies must be used in equation (lO) to reconstruct 
the measured heat inputs, it is consistent also to make use of the wall- to 
bulk-enthalpy difference, i v - i-^, instead of the usual temperature difference, 

T “ in t]ie ^ efini ti°n of the heat-transfer coefficient (eq. (6)). This 
will define a modified heat-transfer coefficient. 


By analogy to equation (6), a modified heat-transfer coefficient can be 
defined by 


h 1 



The relation between h* and h is from the ratio 


ah 


h* _ T w ~ -Pb 
h " i w - i b 

Equation (l2) has the units of reciprocal specific heat. 


( 12 ) 
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The Nusselt number with the modified heat-transfer coefficient is now 
defined as 



where the thermal exchange coefficient 



The modified Du 1 can be used in the Dittus-Boelter equation 


( 13 ) 


(14) 


Nil' = a Re b Prf 
e 


(15) 


where all transport properties must be calculated on the equilibrium basis. 
This correlation still retains the problem of determining the best reference 
temperature for calculating the properties. Besides the reference temperature 
defined by equation (7), reference enthalpies as 


ref 


= i, + m(i 
b w 


V 


(16) 


can be used to estimate a mean reference or film temperature. An analytical 
study on variable property turbulent heat transfer to nonreactive gases showed 
that equation (l6) can be used very effectively with the usual Nu-Re corre- 
lation (ref. 12) . 


Brokaw (ref. 9) suggested the use of mean integral properties ^ which were 
calculated from the expression 


cp = 



cp(T) dT 


(17) 


where cp(T) represents any thermal or transport property. It was found that in- 
tegral mean properties provided the best correlation for the data of Beal and 
Lyerly (ref. 8). 

The data of this report was correlated by the Dittus-Boelter relation with 


13 


I 



( 18 ) 


both standard and modified Nusselt numbers. The rearranged equation 



a Re° 


was the form used to test the data. The value of the Prandtl exponent in equa- 
tion (15) has been taken to be l/3 in equation (18) . 


For correlations with the modified Nusselt number the Reynolds number was 
taken as 



(19) 


where the viscosity was calculated at a specified reference temperature. For a 
Reynolds number defined by equation (l9), the modified Stanton number has the 
particularly simple form 


St’ 


Hu 1 
Re Pr 


or 


St’ 


hi 

G 


( 20 ) 


The St’ of equation (2 0) has a definite advantage as a heat-transfer parameter 
since it contains no terms whose values depend on the choice of the film temper- 
ature. 

From equation (l5). 


St' = a Re b_1 Pr, 


-2/3 


( 21 ) 


Colburns' analogy (ref. 13, p. 424) follows directly from this 


y * St’ Pr 
H e 


2/3 _ 


= a Re 


b-1 


f 

t ' a — 


"H 2 

where the Blausius friction factor expression 

f = 0.079 Re -1 / 4 
was used in the present work (ref. 13, p. 17l) . 


( 22 ) 

(23) 

(24) 
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RESULTS AND DISCUSSION 


The correlated test results of the air data are shown in figure 10. The 
purpose of these runs was to check out the instrumentation and the methods of 
data reduction. 

Plots of wall temperatures for one air run and one run of the NgO^ system 

are shown in figure 11 . This pair of results was chosen because* of similar mass 
rates of flow. It illustrates the profound effect of chemical reaction on heat 
transfer to a gas, for the frozen specific heat of the gaseous N^O 4 system is 

essentially that for air at moderate temperatures, and the system would have 
similar temperature curves to the air run were it not for the presence of the 
reaction N 2 C >4 t- 2N0 2 . In this case, the N 2 C >4 system absorbed about twice the 
heat flux as the air system, and with a significantly lower wall temperature 
distribution and bulk temperature rise. 

A criterion for fully developed heat transfer is a constant wall- to bulk- 
enthalpy difference. This is demonstrated in figure 12 for the same NgO 4 run 

that was used in figure 11. It was obviously impossible to choose, from the 
(T^ - T-^) curve, a point on the tube where fully developed heat transfer was 

established. The enthalpy difference curve indicated that fully developed heat 
transfer was established at the thermocouple position at 6.75 inches on the test 
section. This point corresponded to an x/D = 7.69. The constant property 
analysis of reference 21, for a Reynolds number of 50 000, found that a heated 
length of about 12 tube diameters was required to bring the local Nusselt num- 
ber to within 5 percent of the fully developed value, and about 25 diameters to 
approach within 1 percent of fully established heat transfer. It should also 
be noted that the longitudinal temperature gradients at the entrance to the 
heated tube (fig. ll) were not as severe for the reacting system and affected a 
shorter length of tube than did the air runs. 

A conventional heat-transfer correlation of the experimental data is given 
in figure 13. These are the runs of table II. The abnormally large Nu arises 
from the use of the low values of k f and the anomalies of the conventional 

heat-transfer coefficient defined by equation ( 6 ) . Since figure 12 demonstrated 
that a conventional heat-transfer coefficient would not be constant in the re- 
acting system, it was necessary to arbitrarily choose a thermocouple position 
within the fully developed region from which the necessary temperature data 
could be obtained; in figure 13 all data were computed from temperatures at the 
tube position of 17.63 inches (x/D = 34.6). 

Figure 13 also indicates that the conventional correlation deviates more 
from the analytical result of reference 21 as the experimental wall temperatures 
diminished. This is in part due to the approach of the proper equilibrium 
thermal conductivity k e to the limiting value of k^ at higher temperature 
(fig. 6 ) and in part due to the mean specific heat 
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( 25 ) 



approaching the frozen for large T w as in figure 5. 

Figures 14 to 16 illustrate the variation of equilibrium specific heat 
that occurred in the data for the three levels of wall temperatures. This 
spread of property variables over the three groups of experiments was deliber- 
ately planned, as mentioned in the INTRODUCTION, to give the widest possible 
variation in chemical effects at the wall. 

The conventional correlation of the data with the equilibrium thermal con- 
ductivity k e used in the Nusselt number has been given in figure 17. There 
the results for each run have been given at tube positions of 13.75 and 17.63 
inches (x/D = 25.0 and 34.6, respectively). The contrast between these results 
and those in figure 13 is quite clear: the Nusselt numbers with k g calculated 

at the wall temperature and those with k e calculated at the bulk temperature, 
are essentially mirror images with respect to the region of the two analytical 
curves for constant properties. The choice of wall or bulk temperature made 
little difference in the Nusselt numbers of figure 13, but in figure 17 it had 
a most profound effect. 

The other contrast observed between the two figures was the reversal of 
the positions of the low and high wall temperature Nusselt numbers relative to 
the analytical curve for constant properties. 

The Nu, Re, and Pr functions were also calculated at temperatures corres- 
ponding to reference enthalpy positions of m = 0.3, 0.4, 0.5, and 0.6 in 
equation (16) . The values for m = 0.6 were much closer to the known analyt- 
ical results for constant properties, and they have been plotted in figure 17 
between the analytical curve for Pr = 1.0 (ref. 22), and the curve for 
Pr = 0.7 (ref. 23). The comparison of these data with the analytical curves is 
quite good. 

The basic results from the experimental runs are summarized in table III. 
The heat flux q, the enthalpy difference i w - i fe , and the modified heat- 

transfer coefficient h T , were invariant in the fully developed heat-transfer 
region of the tube, and subsequent correlations have made use of these data for 
this reason. 

The comparison of the modified heat-transfer coefficient h* with the mass 
average velocity for the data of table III, and for the data of references 8 
and 9 is shown in figure 18. Of the data reported in table III, fourteen runs 
were plotted within ±6 percent of the correlation 

h' = 0.00596 G 0 - 884 (26) 
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and of the four remaining runs, the high wall temperature series had deviations 
of 10 and 20 percent, the moderate wall temperature series had one of -22 per- 
cent, and the low wall temperature series had one of 8 percent deviation. 

The data of Beal and Lyerly (ref. 8) and Callaghan and Mason (ref. 5) for 
uniform wall temperature systems have been used for figure 18. Local heat- 
transfer coefficients and data were recalculated to furnish the modified coef- 
ficients used in the figure. The data of Beal and Lyerly show the most devia- 
tion from the correlating curve, and this is probably because their measure- 
ments were taken at the extremely short heating length l/d = 2. 

The success of equation (26) in correlating the is rather sur- 

prising because it does not depend at all on any transport properties such as 
those used in the conventional Dittus-Boelter formulation (eq. (5)). It must 
again be emphasized here that the dissociation data plotted in figure 18 ex- 
hibited the widest possible variation of chemical effects on the properties of 
the system. 

A correlation somewhat similar to equation (26) was independently inves- 
tigated by Simoneau and Hendricks (ref. 24) for data on air, helium, hydrogen, 
and carbon dioxide, with rather large wall- to bulk-temperature differences. 
Their formulation can be written as 



that is, in terms of the standard heat-transfer coefficient (eq. (6)). 

Simoneau and Hendricks found that a different proportionality constant for 
equation (27) was needed with each of the four nonreactive gases. This behavior 
is analogous to the situation between the air and data in figure 18, where 

the form of equation (26) holds also for the air data, but a change in the value 
of the coefficient is necessary to shift the correlating line upward to enclose 
the air data points . 

Correlation of the experimental data by equation (18) has been shown in 
figures 19 to 21 for Pr e , r e , and q evaluated at local bulk and wall temper- 
atures, and as integral mean values between the local wall and bulk tempera- 
tures . 

Of the three correlating curves plotted in the figures, the ones due to 
Deissler (ref. 22) and Sparrow, et al. (ref. 2l) are for constant property 
analyses . 

Callaghan and Mason (ref. 5) found that equation (18) with a = 0.018 and 
b = 0.8 gave a good fit to their data over their limited range of Reynolds 
numbers. It is seen that their equation gave the best correlation for the pres- 
ent data over the extended Reynolds number range for any of the three bases of 
property calculations. 
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The considerable scatter in the wall -based data of figure 20 is due to the 
extreme differences in the wall temperatures between the three sets of data of 
the present experiment. The properties at the wall for the highest T^ runs 
were essentially those of frozen composition, and this was especially true for 
the sections of the tubes where fully developed heat transfer was established. 
Figures 14 to 16 demonstrate that there was far less variation in the bulk 
properties among the three data sets than in the wall properties. This was re- 
flected in the better correlation on a bulk basis as is shown in figure 19. 

The data correlation with the equilibrium properties calculated with equa- 
tion (l7) are shown in figure 21. This correlation is demonstrably superior to 
the previous ones of figures 19 and 20. 

For Reynolds numbers in the range Re = 10^ to 10^ the equation 


Pr' 


Hu* n ft 

-373 = 0.018 Re u ’ y 


(27) 


predicts values for the lu'/Pr^/ 5 term that are about 25 percent lower than 

those using the equation of Sparrow, Hallman, and Siegel (ref. 2l) . This lat- 
ter equation was shown in figure 10 to correlate the air data tests very 
nicely. Thus the air and N^O 4 data discrepancy of figure 18 has been carried 
over into the Nusselt number correlation. A probable explanation of this is 
that these correlations are somewhat inadequate approximations to the basic 
heat-transfer processes that occurred within the reacting system. 


Results of correlating the data by the Colburn analogy of equations (22) 
and (23) are shown in figures 22 to 24. Comparison has been made with the 
analytical data of Deissler, and with the two empirical curves of figures 19 
to 21, all having been transformed to the Colburn factor jjj. With Pr e and 
Re based on bulk or wall temperatures, figures 22 and 23 show considerable 
data scatter in Colburn factors with these bases. Figure 24 shows the correla- 
tion derived with integral mean properties. The agreement with the bottom cor- 
relating curve is fairly good, especially if the four divergent data points 
discussed with figure 18 are discarded. 

The exact Colburn analogy (eq. (23)), is plotted on figures 22 to 24 as 
one-half of the Blausius friction factor (eq. (24)). Deissler r s analysis for 
Pr = 1.0 is most closely approximated by the exact analogy. 


The previous figures show that the experimental Ng O4-NO2 heat -transfer 
data can be correlated with a Dittus-Boelter -Colburn-type relationship and with 
regard to the use of the correct heat-transfer coefficient and thermal prop- 
erties. That the correlation does not conform to the results of constant prop- 
erty analyses is at least in part due to the assumption of the Prandtl number 
entering the relation as a simple power. Sparrow, Halljnan, and Siegel 
(ref. 2l) had earlier concluded on the basis of analytical studies that the 
Prandtl number does not have such a simple power form. Petukhov and Popov 
(ref. 25) came to a similar conclusion after studying a large number of experi- 
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mental data from unreactive gases with property variations. 


CONCLUDING REMARKS 

The following results were obtained in the investigation of local heat- 
transfer coefficients in the turbulent flow of gaseous NgO^-NOg through an 
electrically heated tube: 

1. No consistent data correlation was obtained using the conventional 
Dittus-Boelter heat-transfer equation with frozen properties and with the 
standard heat-transfer coefficient based on the wall- to bulk-temperature 
difference. 

2. Good agreement with constant property analyses was obtained when ex- 
perimental Nusselt numbers were computed from reference equilibrium thermal 
conductivities and the standard heat-transfer coefficient and when reference 
equilibrium Prandtl numbers were used for the system. This correlation was 
very sensitive to the correct reference pointy which was found to be at the 
temperature corresponding to the reference enthalpy ^ref ” 0,6 ^ i w “ ^b^ + ^b * 

3. A simple logarithmic correlation of the experimental modified heat- 
transfer coefficient h ! , which was independent of transport properties, with 
the mass average velocity G was surprisingly good. A similar correlation was 
possible with local air data, but values of h r were about 18 percent higher 
than the dissociation result. 

4. For each experimental run, the local modified heat-transfer coefficient 
h* based on local wall- to bulk-enthalpy differences was constant over the 
fully developed region of the test section. The local standard heat-transfer 
coefficient h based on local wall- to bulk-temperature differences, continu- 
ously decreased along the tube in the flow direction. Conversely, by using the 
constant value h T as the definition of fully developed heat transfer, the 
thermal entrance length for the dissociating system was found to be about half 
that predicted from analysis for constant property fluids. 

5. Local Colburn factors Nu 1 /Prl/3 

were computed from the experimental 

modified heat-transfer coefficients h 1 and properties evaluated at three ref- 
erence points: wall temperature, bulk temperature, and the mean temperature 

integral of the properties. Of these, the integral mean properties provided 
the most consistent data correlation, which agreed well with the expression 

jj-7/pJl/3 = 0.018 Re 0 * 8 . 

6. The Colburn analogy was shown to hold quite well for Deissler’s con- 
stant property analysis with Prandtl number of unity, but the correlation of 
the experimental NgO^-NOg data, with integral mean properties, was quite good 

with the expression St* Pr 2 / 3 = 0.018/Re 0 * 2 which predicts values of the 
Colburn factor about 20 percent less than those of Deissler’s constant 

property analysis. 
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7. The Dittus -Boelter -Colburn heat-transfer correlation of the present 
experimental data, as well as of those of other investigations using N 2 0 4 -N0 2 
gas, was consistently lower in magnitude than the same correlation for constant 
property analytical data. 


Lewis Research Center, 

Rational Aeronautics and Space 
Cleveland, Ohio, November 


Administration, 
2, 1965. 
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APPENDIX 


PROPERTY CALCULATIONS 


Degree of Dissociation 


Determination of the degree of dissociation in the gas at any point in the 
test section was necessary for all further calculations of densities , and 
thermodynamic and transport properties. The basis of the calculation is the 
equilibrium constant expressed in terms of the partial pressures of the react- 
ants and products 2 


*p- 


■NOc 


•NoO 


(Al) 


2 4 


For a Dalton r s Law system the partial pressures are products of the total 
system pressure P and the molar fractions of the constituents 


■NO, 


= y. 


NO, 


and 




3 = y P 

W 2 °4 N 2°4 


(A2) 


where + y = 1 by definition. For an initial 1 mole of N o 0„, the 

NOr- ivr_ n . 17 *-> 4 = 


n 2 o 4 


equilibrium composition has, for a mole fraction of NgO^ decomposed, 
(l - a) moles of NgO^ remaining and 2a moles of NO 2 formed from the 


n 2 o 4 . 

hence. 


There are then (l - a) + 2a = (l + a) moles in the equilibrium mixture; 




2a 


and 


NOc 


% 2 ° 4 


1 + a 


1 - a 
1 + a 


(A3) 


From the aforementioned expressions the equilibrium constant in terms of the 
degree of dissociation is 




p 

4a P 



(A4) 


In terms of the equilibrium constant and the system pressure, the degree of 
dissociation is 
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I 



a 


(AS) 



Equation (A 5) was the basis for all calculations of the degree of dissociation. 


At the standard temperature 298.16° K (25° C) the equilibrium constant can 
be calculated from the standard state Gibbs 1 free energy, since (ref. 26, 
pp. 334, 349) 


AE° 

- = R In Kp (A6) 

The necessary thermodynamic data is given in reference 31. 

Extension of the equilibrium constant calculation to temperatures above 
298.16° K was accomplished by integration of the van’t Hoff isobar (ref. 26, 
p. 338) . 


/d In Kp\ ^ £jjO 

V /p RT 2 

The standard enthalpy of reaction AH° at 298.16° K was obtained from refer- 
ence 27. Variation in the AH° at higher temperatures, due to molar specific 
heat changes, was calculated by the method outlined in the section Enthalpies 
of this appendix. 


(A7) 


The empirical equation 


logpQ Kp — 9 • 218 


3006 

T 


(A8) 


was found to give an excellent fit to the equilibrium constant data up to 
490° K. At higher temperatures the agreement was somewhat less successful, but 
the gas was also almost completely dissociated at these temperatures and for 
the pressures encountered in the experiments. 

The equilibrium constant increased rapidly with increasing temperature 
from a value of 0.157 at 300° K to 1180 at 490° K. Equation (A5) shows that 
large values of IC result in values of a approaching unity. 


Mixture Densities 

An ideal mixture of perfect gases obeys the perfect gas law 

PM 
P RT 


(A9) 
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where 


M 


is the average molecular weight of the mixture, that is. 


for the gas system, 

for the mixture 


” = + y ^2°4 

Since « N2 o t - 2 « H 0 2 - 



(A10) 


equa t i ons ( A3 ) and ( A10 ) g i ve 


n 2 o 4 


1 + a 


(All) 


From equations (A9) and (All), the mixture density is 

p % 2 0 4 

p RT(1 + a) 


(A12) 


Enthalpies 


The ideal gas enthalpies of a nonlinear polyatomic molecule are calculated 

from 


H 


o 



4R + R 




(A13) 


where the Xj = 1.4384 ci)j/T and the summation is over all observed wave num- 
bers (ref. 28). This reference gives the numbers for the N0 2 and molecules 

as listed in table IV. 
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energy for the j « 7 line,, since that wave number is essentially zero in 
value. From this, for Xy ^ 0, 

X. 

lim -^-s2 = 1 

X. ^ o e J “ 1 

J 


The enthalpy of reaction at a temperature T is the difference between 
the summation over the product enthalpies and the summation over the reactant 
enthalpies (ref. 29) . For the system ^ 2M> 2 the molar reaction enthalpy 

is 


AH° = 2H° 
NCL 


- H 


N 2°4 


(A16) 


Expressing the right-hand side of equation (A16) in terms of the respective 
relative enthalpies (ref. 29) gives the molar enthalpy of reaction as 


0 


AH° « 2 |(H° - H°) + AH° 


- I 

jCjNOg 


(H° - H°) + AH° 
v O' f 




= 2(H° - H°) - (H° - H°) + 2(AH° ) - (AH° ) (A17) 

0 U0 2 0 N 2 °4 f '° N0 2 ?>0 N 2 0 4 

The molar heats of formation at 0° K are given in reference 27 as the following 


for WOg : 


and for 


AH° = 8.682 kcal/mole 
f,0 


AH° n = 4.489 kcal/mole 

I ^ u 


The (H° - H°) terms in equation (A17) are calculated with equations (A14) and 
(A15) . Enthalpies of the gas mixture on a -unit mass basis were calculated from 


mix ^m 2 x m 2 


+ Jl 


N 2°4 lN 2°4 


(A18) 


where in terms of the degree of dissociation of V-,0^ the mass fractions of 
the species are 





w 2 o 4 


1 - a 


(A19) 


The specific enthalpies of the individual species (i.e., enthalpies on a unit 
mass basis) are readily obtained from equations (A14), (A15), and (A17) as 
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and 


H- 




N0 2 M- 


'NO; 


Mi 


NO; 


( 


(H° - 


H°) + AH° 
O' f 




no„ 


H- 


TJ 2 °4 


NO 
2 4 


M 


N 2°4 


M- 


N 2°4 


(h° - 


H°) + AH° 1 

0 ^oJ H2 o 4 


(A20) 


(A21) 


Equation (A18) was the basis for enthalpy calculations with the experiment 
together with equations (A19) to (A2l) . 


Specific Heats 

The specific heat on a mass basis of the reacting gas mixture is derived 
from the specific enthalpy (eq. (A18)) in the usual manner 


S . 

p,e \ ST / 


(A22) 


'di 


NO, 


'hi 


a | 


j^>4 

St 


(l - a) + i 


N0 2 


i „ 2 0 4 (I) 

2 4 p 


(A23) 


The terms on the right-hand side of equation (A23) are evaluated from equations 
(A19) to (A2l) . This results in 


(I) 


^ -N # . ^ AH T> /Xy> 

C p,e " ^p.NOg + (1 - a)C P,W 2 0 4 + ^ 


%,e ~ °p,f + °p,r 


(A24) 

(A25) 


where the "frozen" (i.e., constant composition) specific heat is 
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I 



(A26) 


G p,f " aC p,N0 2 


+ (1 - a)C p ^ 2 o 4 


and the reactive specific heat is 


; _ AH° /cia\ 

' p ' r " % 2 0 4 W P 


(A27 ) 


The gradient of the degree of dissociation with temperature at constant 
pressure uses equations (A5) and (A7) for the derivation. It can he shown from 
equations (A4) and (A7) that 


5a 

5T 


54 2\ AH 

-n 1 


RT 


(A28) 


so that the reactive component of specific heat becomes 


P, r 2 M n 2 0 4 


ccA - a 2 ) 

i\r„n, \ / 


(A29) 


Equations (A24) and (A28) were those used to compute the experimental specific 
heats . 


Thermal Conductivity- 

Experimental thermal conductivities of gaseous N^O^ have been measured 
by Coffin and O'Neal (ref. 3) . These values were shown to agree with analyt- 
ical data calculated with the method of Butler and Brokaw (ref. 2). These 
authors showed that the thermal conductivity of a gaseous system in chemical 
equilibrium could be expressed as the sum of a frozen conductivity and a re- 
active conductivity ^ the latter representing the diffusional transport of the 
enthalpy of reaction 


k e = k f + k r 


Reference 2 showed that the reaction term 


(A30) 


= 




(A31) 


The term in brackets is independent of pressure and its values between 290° K 
and 490° K have been tabulated by Coffin and O’Neal (ref. 3). 

The frozen (i.e., constant composition) contribution to thermal conduc- 
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tivity is expressed as a sum of two terras 


k f = k» + k” (A32) 

where k* accounts for the collisional transport of translational kinetic en- 
ergy of the species in the mixture and k" accounts for the diffusional trans- 
port of internal energy. The forms of these two terras are somewhat involved, 
and references 30 and 31 should be consulted for details. Reference 30 is 
especially recommended, although it is based on the treatise of Hirschf elder , 
et al. (ref. 3l) . The frozen contribution to thermal conductivities needed for 
this experiment were calculated from the equations of reference 8, and the 
force constant data for the NOg and NgO^ molecules were taken from table I 
of that reference. 

The reactive contribution was calculated with equation (A3l) . Values of 
the term (®p/rt) (AH o2 /rT 2 ) in equation (A3l) were interpolated from the data 
of Coffin and O’Neal (ref. 3). In order to extrapolate their data beyond 
490° K the following empirical expression was used: 

•0PAH° 2 2.756 . . 

RT 2 ~ 1.201 

RT T 


Viscosities 

The viscosity of the gaseous NgO^ system was calculated with equation (13) 
of reference 30. The form of this equation is similar to that for computing 
the k M term in the frozen thermal conductivity. As with the thermal conduc- 
tivity calculations, the force constants of the N0 2 and N 2 0 4 molecules were 
obtained from table I of reference 30. 

Three sets of data are recorded on the viscosity plot of figure 7. The 
data of reference 7 is most certainly in error. The more recent data of Beer 
(ref. 32) shows the most discrepancy at the lower temperatures; the data of 
Petker and Mason (ref. 33) collected with a rolling-ball viscometer, parallels 
the curve of Brokaw (ref. 18), but with a discrepancy of 15 percent at 450° K. 
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TABLE I. - TEST CONDITIONS FROM VARIOUS DATA SOURCES 


Source 

Type 

Boundary 

condition 

(a) 

Tube length 
to diameter 
ratio, 

l/d 

Heat flux 
Btu/(hr) (sq. ft) 

Entrance 

Reynolds 

number, 

Re 

Maximum 
local wall 
temperature, 
°K 

Inlet 

pressure, 

atm 

Reference 4 

Experimental 

UWT 

81.2 

4.1 to 6.5X10 3 

12.0 to 20.0X10 3 

375 to 387 

1 . 

0 

Reference 8 

Experimental 

UHF 1 

2 

2.6 to 12.5 

30.7 to 111 

315 to 374 



Reference 9 

Analytical 

UHE 

2 

2.6 to 12.5 

30.7 to 111 

315 to 374 



Reference 5 

Experimental 

UWT 

38.0 

1.7 to 5.9 

11.3 to 19.9 

337 to 442 

* 

f 

Reference 6 

Experimental 

UWT 

88.0 

.3 to 8.4 

6.7 to 22.0 

311 to 477 



Reference 7 

Experimental 

UWT 

16.0 

.3 to 13.3 

26.8 to 209 

321 to 461 

0.9 to 1.0 

Reference 10 

Analytical 

UHF 



10.0 to 200 

300 to 370 

1 . 

o 

Reference 11 

Experimental 

UHF 

184.0 

.6 to 27.0 

5.6 to 68.2 

315 to 450 

0.9 to 1.7 

Present 

Experimental 

UHF 

39.6 

6.2 to 44.2 

18.9 to 135 

359 to 572 

1.6 to 2.1 

investiga- 









tion 










cl 

UWT = uniform vail temperature. 
UHF = uniform heat flux. 



TABLE II. - EXPERIMENTAL CONDITION'S FOR N 2 0 4 SYSTEM 


Run 

Pressure, 

Flow rate. 

Heat 

Bulk 

Maximum 

Heat 


p, 

w. 

input , 

temperature 

rail 

balance 


atm 

lb/hr 

Btu/hr 

rise, 

V 

°K 

temperature, 

V 

°K 

error, 

percent 

1 

1.9 

58.90 

3995 

344 to 365 

562 

11.1 

2 

1.9 

53.20 

4165 

330 to 353 

562 

4.3 

3 

1.8 

44.40 

3960 

314 to 343 

562 

2.6 

4 

1.8 

34.00 

3228 

312 to 343 

568 

3.9 

5 

1.8 

23.65 

2422 

312 to 345 

572 

5.2 

6 

1.7 

14.80 

1663 

309 to 346 

566 

6.4 

7 

1.9 

89.40 

2808 

318 to 329 

359 

14.9 

8 

1.7 

64.10 

2388 

320 to 331 

362 

11.4 

9 

1.8 

66.20 

2148 

320 to 330 

359 

8.8 

10 

1.8 

67.40 

2198 

319 to 330 

360 

6.4 

11 

1.7 

47.80 

1846 

316 to 329 

360 

4.5 

12 

1.6 

31.30 

1403 

315 to 330 

364 

.2 

13 

1.6 

17.93 

874 

313 to 329 

361 

3.0 

14 

2.0 

115.8 

6240 

326 to 340 

447 

12.7 

15 

2.1 

89.10 

4750 

325 to 342 

440 

1.5 

16 

2.1 

75.50 

4243 

325 to 342 

444 

1.0 

17 

2.1 

66.80 

3836 

325 to 343 

445 

2.8 

18 

2.1 

46.20 

2892 

321 to 342 

447 

3.4 



TABLE III. - SCME EXPER IMEHTAL RESULTS FOR EgOq SYSTEM 


Run 

Mass 

Uniform 

Fully 

Modified 

Experimental 

Inlet 


average 

heat flux, 

developed 

heat- 

modified 

Reynolds 


velocity, 


enthalpy 

transfer 

Stanton 

number, 


G, 

, 

Btu/(hr) (sq ft) 

difference, 

coefficient. 

number, 

Re 


lh/(sq ft) (hr) 

\ " V 

Btu/lb 

hS 

h* 




lb/ (sq ft) (hr) 

G 


1 

66 200 

28 520 

155.4 

178 

2.69X10" 3 

61 . 3X1 0 3 

2 

59 760 

29 530 

195.7 

151 

2.53 

60.4 

3 

49 900 

28 080 

229.0 

123 

2.46 

55.3 

4 

38 200 

22 880 

232.4 

98.3 

2.58 

42.9 

5 

26 580 

17 170 

229.3 

CD 

,, sH 
t — 

2.81 

29.9 

6 

16 630 

11 790 

227.3 

51.9 

3.12 

18.9 

7 

100 400 

19 910 

104.0 

192 

1.91 

108.7 

8 

72 000 

16 940 

98.96 

171 

2.38 

77.1 

9 

74 400 

15 230 

93.91 

162 

2.18 

79.8 

10 

75 750 

15 690 

95.59 

163 

2.15 

81.6 

11 

53 720 

13 080 

100.9 

130 

2.42 

58.6 

12 

35 180 

9 945 

107.3 

93.4 

2.65 

38.6 

13 

20 140 

6 196 

102.3 

60.6 

3.01 

22.4 

14 

130 100 

44 250 

175.1 

253 

1.94 

135.1 

15 

100 100 

33 680 

174.9 

192 

1.92 

104.6 

16 

84 800 

30 100 

176.3 

171 

2.02 

88.8 

17 

75 050 

27 200 

175.8 

155 

2.07 

78.6 

18 

51 920 

20 520 

182.1 

113 

2.18 

55.4 











Figure 1. - Schematic of test loop. 







(a) Test section. 


/-Pressure tap 



(b) Upstream electrode flange. 

Figure 2. - Test section detail and thermocouple location. (All units are inches. ) 



Temperature, °K 

Figure 4. - Equilibrium and frozen specific heats as function of tem- 
perature for N 0 O 4 system. Absolute total pressure, 2 atmospheres. 




Thermal conductivity, k, Btu/(hrHftK°R) 



Temperature, °K 


Figure 5 . - Equilibrium and frozen thermal conductivities as function 
of temperature for N2O4 system. Absolute total pressure, 2 at- 
mospheres. 



Temperature, °K 


Figure 6. - Specific enthalpy as function of temperature for N2O4 
system. Absolute total pressure, 2 atmospheres. 



Figure 7. - Dynamic viscosity of - 2N0 2 system as function of 
temperature for absolute total pressures of 1 and 2 atmospheres. 



Figure 8. - Thermal exchange coefficient as function of temperature 
for N 2 0 4 system. Absolute total pressure, 2 atmospheres. 
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Figure 9. - Equilibrium and frozen property Prandtl numbers as func- 
tion of temperature for ^4 system. Absolute total pressures, 1 and 
4 atmospheres. 



Figure 10. - Conventional heat-transfer correlation for air data. Properties 
evaluated at reference temperature T ref = T b + 0. 5<T W - T b ). 
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Figure 12. - Variation in wall- to bulk-enthalpy and temperature differences along tube for N 2 O 4 system. 
Run 4. 
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Figure 13. - Conventional heat-transfer correlation for N 2 O 4 system 
with frozen properties. Prandtl number Prf = 0. 708 to 0. 713. 


Enthalpy difference, i w - i b , Btu/lb 






Equilibrium specific heat, C D e , Btu/(lb)(°R) 



Figure 16. - Wall to bulk variation in equilibrium specific heat for high 
wall temperature runs in N 2 O 4 system. Runs 1 to 6 . 



Nusselt number, Nu = hD/(k p ) 



Figure 17. - Conventional heat-transfer correlation of N 2 O 4 data with equilibrium properties. 
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4x10-’ 


Figure 18. - Basic correlation of several N 2 O 4 heat-transfer experi- 
ments. 



Figure 19. - Data correlation for N 2 O 4 system with bulk properties. 
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0 oo 



Figure 20. - Data correlation for system with wall properties. 



Figure 21. - Data correlation for N 0 O 4 system with integral mean properties. 



Colburn factor, (ju) = St'<Pr p ) 



Reynolds number, Re b = GD 

Figure 22. - Correlation of N 2 0 4 data by Colburn analogy. Properties evalu- 
ated at bulk temperatures. 



Figure 23. - Correlation of N 2 O 4 data by Colburn analogy. Properties evaluated at wall 
temperatures. 
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Figure 24. - Correlation of N^data by Colburn analogy. Properties evalu- 
ated at integral mean. 


NASA- Langley, 1966 E-2893 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 
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